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Abstract. Magnetoelectrics often possess ions located in noncentrosymmetric 
surroundings. Based on this fact we suggest a microscopic model of magnetoelectric 
interaction and show that the spin-orbit coupling leads to spin-dependent electric 
dipolc moments of the electron orbitals of these ions, which results in non-vanishing 
polarization for certain spin configurations. The approach accounts for the macroscopic 
symmetry of the unit cell and is valid both for commensurate and complex 
incommensurate magnetic structures. The model is illustrated by the examples of 
M11WO4, MnPS3 and LiNiP04. Application to other magnetoelectrics is discussed. 
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1. Introduction 

Pierre Curie pQ was the first to predict the interplay between magnetism and 
ferroelectricity. The macroscopic symmetry consideration of the magnetoelectric (ME) 
effect was given by Dzyaloshinskii only in 1959 [2], whereas experimentally it was 
discovered by Astrov in Cr 2 03 in 1960 |3J. In the last decade whole new classes of 
magnetoelectric materials were discovered and experimentally studied (for a recent 
review see, for example, jl]). 

Despite a long history of ME effect its microscopic origin is still a subject of 
debate. The fact that ferroelectricity commonly occurs in incommensurately modulated 
magnetically ordered phases predominates the assumption and consideration of various 
complex types of magnetic order such as screw, cycloidal, helix and others in most 
of the microscopic models of magnetoelectricity [U [5]. From the macroscopic crystal 
symmetry point of view the close connection between the appearance of modulated and 
ferroelectric phases in magnetoelectrics was recently pointed out [6]. Whereas electric 
polarization indeed often occurs concurrently with complex modulated spin structures 
it can also be induced by commensurate magnetic order (such is the case, for example, 
in some rare-earth manganates RM^Os [7]). 

Currently two models of ME coupling are widely accepted in literature. In the 
model by Sergienko et al. {8] electric polarization P ~ [Si x Si+i] is induced by the 
Dzyaloshinskii-Moriya (DM) interaction between two magnetic ions % and % + 1 with 
superexchange. In the spin current model by Katsura et al. [9] the electric polarization 
is due to the spin supercurrent P ~ x [Si x Sj) with being the unit vector connecting 
the sites % and j. Interpretation of the experimentally observed polarization using these 
models meets difficulties though and was a subject of critique pllj EH]. It was argued [H] 
that the DM coupling [8] is 2 orders of magnitude weaker than what is needed to explain 
experimental situation in magnetoelectric manganites. On the other hand in the spin 
current model [9] the authors also overestimated up to 2 orders of magnitude the value 
of the magnetically induced polarization [TT] . Several other mechanisms of the ME effect 
were proposed [TTJ [12] but as argued by the authors themselves they are insufficient to 
explain magnetoelectricity in some Cu 2+ magnetoelectrics. 

From our point of view most of the microscopic models proposed so far do not take 
into account the macroscopic crystal symmetry mainly focusing on three-site clusters 
(two metal ions and oxygen) and considering various spiral magnetic structures. The 
authors of the DM model of ME effect [8], for example, apply it to the case of rare- 
earth manganites RMnOs, whereas it can be shown that the macroscopic polarization 
vanishes when one considers the orthorhombic symmetry of the unit cell. Indeed, the 
authors take TbMnOs as an example using the spiral magnetic structure available from 
neutron diffraction data. They consider the polarization induced in the x and y chains 
of Mn-O-Mn, but erroneously take into account only the z = plane. Indeed, when one 
considers also the z = 1/2 plane (obtained from the first one by cr 2 (00|)) the macroscopic 
polarization cancels out. 
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Figure 1. (a) The monoclinic unit cell of M11WO4. Mn 2+ , W 6+ and 2 ~ ions are 
shown by black, dark gray and light gray circles, respectively, (b) The Mn 2+ cluster 
showing distorted oxygen octahedron. The rotational C2 symmetry axis is indicated. 

The aim of any ME model is to find the mechanism of inversion symmetry 
breaking by magnetic order. At the same time some authors explicitly start with 
a centrosymmetric cluster such as, for example, in the spin current [HI [13] or the 
MeO n cluster model [TT] . This dismisses the fact that in many magnetoelectrics 
some types of ions are located in noncentrosymmetric surroundings already in the 
paramagnetic phase. In this work using MnWC>4 as an example we suggest a microscopic 
magnetoelectric coupling model taking into account macroscopic symmetry of the 
unit cell and noncentrosymmetric surrounding of the Mn 2+ ions. We then apply our 
microscopic model to estimate the linear magnetoelectric coefficients in M11PS3 and 
LiNiPCU and discuss its application to other magnetoelectrics. 

2. Microscopic model 

Wolframite MnWC>4 possesses a monoclinic structure at room temperature (figure 
described by the space group P2/c (C| ft ). On lowering the temperature it undergoes 
a sequence of magnetic phase transitions at 13.5 K (Tjv), 12.7 K (T 2 ) and 7.6 K 
(Ti), leading to the appearance of magnetically ordered states AF3, AF2 and 
AF1 [14], respectively. The low temperature phase AF1 is characterized by the wave 
vector k = (1/4; 1/2; 1/2), whereas the incommensurate phases AF2 and AF3 by 
(—0.214; 1/2; 0.457). Electric polarization in MnWC>4 appears in the AF2 phase along 
the crystal b axis [15]. The phenomenological model of phase transitions in wolframite 
was suggested earlier based on the assumption that the magnetic order is driven by 
the instability in the (1/4; 1/2; 1/2) point of the Brillouin zone |16J . However, for 
our purpose of building a microscopic model of the ME effect we first start with a 
hypothetical magnetic order with k = in MnWC>4 (i. e. without multiplication of the 
unit cell) and then consider the real magnetic structure. 

MnWC>4 contains two magnetic Mn 2+ ions Mni and Mn 2 in the unit cell located 
at positions (0.5; 0.6853; 0.25) and (0.5; 0.3147; 0.75), respectively. In the following we 
define the orthogonal x, y and z axes parallel to the a axis, parallel to the b axis 
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Table 1. IR's of the P2/c space group corresponding to k — 0. The last column lists 
the components of F, A and P according to the IR's upon which they transform. Note 
that F and A are odd under time inversion. 



IR 


C 2y (OOi) 


I (000) 


Order parameters 


GM 1+ 


1 


1 




GM 1 - 


1 


-1 


Py, Ay 


GM 2 + 


-1 


1 


F F 
± X-) ± Z 


GM 2 - 


-1 


-1 


Px i Pz i A x , A z 



and perpendicular to both the a and b axes of the monoclinic cell, respectively. For 
the description of the magnetic order with k = we can introduce ferromagnetic and 
antiferromagnetic order parameters F = S\ + S 2 and A = Si — S%, respectively, where Si 
and S2 are the magnetic moments of Mnx and Mn 2 . Table [1] summarizes the irreducible 
representations (IR) possessed by the P2/c space group in the center of the Brillouin 
zone and shows how the components of F, A and electric polarization P transform 
according to the symmetry operations of the group. As evident from the table the 
following magnetoelectric interactions are allowed by the macroscopic symmetry 
PAP PAP PAP PAP 

1 Il-^x*- yi 1 [i^iy- 1 - xi 1 n^iy- 1 - zi 1 fi-^z- 1 - yi (\\ 

P A F P A F P A F 

1 y-^a- 1 - ay 1 y-^x- 1 - zi 1 y-^-z 1 xi 

where fj, = x,z and a = x,y, z. Thus any combination of ferromagnetic and 
antiferromagnetic ordering produces electric polarization. The polarization also arises 
in case when only one of the Mn 2+ moments orders (e. g. when Si 7^ and S2 = 0, 
which corresponds to F — A). Such collinear commensurate magnetic ordering gives 
no polarization according to any of the microscopic models proposed so far, since they 
are usually concentrated on noncollinear modulated magnetic structures. Nevertheless, 
from the crystal symmetry point of view such magnetic order breaks the crystallographic 
equivalency of Mni and Mn2 atoms, which are connected by inversion (i. e. breaks 
inversion symmetry) and electric polarization arises. 

In order to build a valid microscopic model of ME interaction we proceed with the 
following consideration. First we note that the Mn 2+ ions in MnWC>4 are located in 
noncentrosymmetric positions. This can be seen from the fact that I (000) interchanges 
Mni and Mn2 or by direct examination of oxygen positions around Mn 2+ ions [14J. 
Since C 2 y (00|) transforms each Mn 2+ ion into itself their local crystal symmetry is C2 
as depicted in figure [Tb. Thus the manganese ions are located in polar surroundings 
and all their electron's orbitals have electric dipole moment along the y axis. Every 
monoclinic cell has two Mn 2+ ions with opposite dipole moments, thus, preserving zero 
dipole moment of the unit cell as schematically shown in figure |2K 

Figure [2b shows an ion in polar crystal environment (local symmetry C2). All 
electron orbitals possess, therefore, an electric dipole moment. If one takes the spin 
degree of freedom into account then the spin-orbit coupling alters the electric dipole 
moment. When the spin directs away from the two-fold rotational axis it breaks the 



Magnetoelectric effect due to local noncentro symmetry 



5 




Figure 2. (a) Two Mn 2+ ions in the monoclinic unit cell with opposing electric 
dipole moments (gray arrows), (b) An ion (black circle) in C2 crystal field (indicated) 
possessing dipole moment do (gray arrow), (c) Same as (b) but taking the spin degree 
of freedom into account. A spin S (white arrow) deviating from the C2 axis breaks the 
rotational symmetry and changes the electric dipole moment to d deflecting it from 
the initial symmetry axis, (d) Two Mn 2+ ions in the monoclinic unit cell with spin 
dependent electric dipole moments d\ and c?2 resulting in nonzero total dipole moment 
of the unit cell d\ + di . 



rotational symmetry resulting in modification of the dipole moment as schematically 
shown in figure [2b. We now proceed from this simple physics consideration to 
semiquantitative quantum mechanical treatment. 

The manganese ions in M11WO4 are located in distorted oxygen octahedra as shown 
in figure [lb. Therefore we start with a Mn 2+ ion in octahedral crystal field and for 
simplicity consider only the e g orbitals. Inclusion of t 2g levels is straightforward. Thus 
as a zeroth order perturbation we have 

where 7 = z 2 or x 2 — y 2 , H is the hamiltonian including octahedral splitting field and 
Ed is the e g energy level. Next, we include as perturbation the monoclinic crystal field 
of C2 symmetry assuming the polar axis along the z axis 



Vcf = c z z + c xy xy + c' xy { 
and spin-orbit coupling 

V so = -KL-S). 



x 



y 



+ d z z 2 



(2) 



(3) 
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Here c z , c xy , c' xy and c' z are coefficients, L is the angular momentum operator, S is the 
spin and A is the spin-orbit coupling constant. Thus, the perturbed hamiltonian has the 
form H = H + V, with V = V CF + V SO - 

The perturbation V mixes the unperturbed 3d e g states with other states and for 
simplicity it is sufficient to consider only the 4p states H \p a ) = E p \p a ) with the energy 
E p , a — x, y, z. The two-fold degeneracy is removed in the first order and one can write 
for one of the perturbed eigenvectors 

h/>)H0)+$>aba>, (4) 
a 

where A a are coefficients and |0) is an eigenvector from the subspace spanned by |<i 7 ). 
The electric dipole moment is given then by 

d= (ip\e.r\ip) = ^2 A a (0\er\j) a ) + c.c. 

a 

Assuming |0) = q\\d z 2) + q2\d x ' 1 -y 2 ) we obtain in the first order of perturbation the 
z-component of the electric dipole moment induced by the local crystal field 

2qlc z et 2 zz2 

«0z — ^ ) 

where A = Ed — E p and t an = (p a \a\d^). Due to local C2 crystal field symmetry 
do x = d 0y = 0. Performing the perturbation up to the third order we get the spin- 
independent first order contribution and the spin-dependent part of the dipole moment 

as 

/ A\ 2 
d x = ~Qx I ^ I S X S Z , 

dy = —Qy \ ) S y S z , (5) 



A x 2 



where 



d z = d 0z + d 0z l-j (s 2 x + s;), 



2q 1 c z et z , z 2 (q 1 t a , z 2 + q 2 t a ^ x 2 



A 

Thus, in addition to the crystal field induced electric dipole moment the spin-orbit 
coupling gives rise to spin-dependent contribution. 

We now use the obtained results for the above case of magnetic ordering in MnWC>4 
with k = 0. Each unit cell has two Mn 2+ ions in opposite polar surroundings with polar 
axes along y. Therefore, using ([5]) for both of them and performing the proper cyclic 
permutation of indices x, y and z we get the electric polarization 

P* = ~Qx^ [^j (A y F x + A x F y ), 

Py = d 0y ^(^J (A X F X + A Z F Z ), (6) 

P Z = -QZ^(^) (AyF z + A Z Fy), 
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Table 2. Components along x of Mni and Mn2 spins induced by {i]i x ,^ix) and 
(r]2x, £,2x) in the reference unit cell (0 • a±) and in the unit cell displaced by one period 
of translation along a (1 • a%). The last two lines show the spin structures in each cell 
rewritten in terms of F and A. 
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■ V2x ~ 
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??2x 4 


- £lx 4 &x 




Vix 4 


- V2x 4 


€lx 


— &x 
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where v is the unit cell volume. Therefore, we obtain the same polarization as the one 
implied by the ME interactions ([I]). Figure EH schematically illustrates relations (jHJ) 
showing Mn 2+ ions with spins Si and 5*2 directed in such a way that the spin-dependent 
electric dipole moments d\ and d 2 result in nonzero total dipole moment of the unit cell 
di + d 2 . 

Other ME interactions present in (CQ) and absent in (jHJ) can be obtained by 
performing the quantum perturbations to higher orders. It has to be noted, that when 
deriving the spin-dependent dipole moments (jSJ) we assumed the polar distortion of 
the octahedral crystal field to be directed along the z axis (i.e. with two oxygens on 
the symmetry axis), which is not the case in MnW04, but this does not change our 
semiquantitative model. 

3. Magnetic ordering in MnW0 4 with k = (1/4; 1/2; 1/2) 

Equations (JB} give the electric polarization for magnetic structures with k — 0. In 
MnWC>4 the magnetic phase transitions can be described by the order parameters with 
k = (1/4; 1/2; 1/2) [16] . In this point of the Brillouin zone the space group P2/c 
possesses two two-dimensional IR's Gi and G 2 - Using the magnetic representation 
analysis conducted earlier [16J, for every direction a = x,y,z we introduce two order 
parameters (r)i a , £i a ) and (rj2a,^2a) transforming according to G\ and G7 2 , respectively. 
The spin components along x induced by these order parameters are given in table EJ 
For every unit cell one can rewrite the spin components in terms of F and A as shown 
in the table. Analogous analysis can be conducted for the y and z spin components. 

The magnetic unit cell in MnWC>4 is 16 times the crystallographic one. We now 
sum up A X F X over the magnetic cell ^ A X F X = 2(rji x T]2x J r^ix^2x) t° obtain the electric 
polarization using (EJ) 

Py = d 0y ^ ; (rjlxVZx + flaSte)- (7) 

Equation (j7|) gives P y in accordance with the polarization that can be derived from 
the ME interaction P y (r]i x r]2 X + £,i x £,2x) obtained from the macroscopic symmetry 
analysis [T6] . In agreement with the experiment [TJ)] and phenomenological model [TH] P y 
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arises in the AF2 phase when both G\ and G% condense. Long-wavelength modulation 
does not lead to cancellation of (J7j) as shown earlier [16]. Similar to ([7]) other 
contributions to polarization can be obtained using remaining ME interactions (fl]). 

The numerical value of P y from ([7]) for MnW0 4 can be estimated as follows. We 
use the real crystal data and oxygen positions from [T4] to perform the crystal field 
expansion (T5]) and obtain c z ~ 4.9 ■ 10~ 9 N and v ~ 138 A 3 . For the matrix elements 
t an we use the hydrogen-like orbitals and obtain t ZjZ 2 w 0.67a /Z, where a is the 
Bohr radius and Z is the charge of the nucleus and core electrons in units of e. Using 
Z f» 5, A ~ 0.05 eV, A 1 eV and noting that qi ~ 1, (r]i x r] 2x + ^i^x) ~ 1 we 
obtain P y ~ 17 /iC/m 2 in good agreement with the experimental value of the order of 
50 /iC/m 2 [15], HZ] - In this numerical estimate, however, we considered only one e g level 
whereas other 3d electrons give comparable contributions to P y . 

4. Discussion 

Thus, we have built a microscopic model of ME interactions in magnetoelectrics using 
MnW0 4 as an example. Starting with the hypothetical magnetic order with k = 
we determined the ME interactions ([I]). Noting that Mn 2+ ions in wolframite are 
located in noncentrosymmetric polar surroundings we suggested a microscopic model 
of magnetoelectricity. The orbitals of 3d electrons of Mn 2+ ions possess electric 
dipole moments due to the crystal field influence, which gives rise to additional 
contributions (such as the considered 4p states) to their wave functions (TjJ. The spin- 
orbit coupling ([3]) induces spin-dependent electric dipole moments ([5]) since the angular 
momentum operator L mixes different \p a ) states. In this part our approach combining 
phenomenological and microscopic models resembles that suggested for the description of 
weak ferromagnetism by Dzyaloshinskii [18] and Moriya [19] . Our approach differs from 
those of Sergienko et al. [8] and Katsura et al. [9] who obtain the electric polarization 
as either P ~ [Si x 5i+i] or P ~ x [Si x Sj], which is essentially a result of 
interacting spins located at two different ions. On the contrary, in our model due to 
local noncentrosymmetry each magnetic ion has spin-dependent electric dipole moment. 
We show that for certain spin configurations the sum over all magnetic ions of these 
spin-dependent electric dipole moments gives rise to macroscopic polarization. 

We then apply this approach to describe real magnetic structures in MnW04 with 
k = (1/4; 1/2; 1/2) by summing up local contributions to polarization in every unit cell 
of the magnetic cell. Our approach naturally accounts for the macroscopic symmetry 
of the unit cell and is valid in both cases of collinear (hypothetical) magnetic ordering 
with k = and complex long-wavelength modulated magnetic structure observed in 
MnW0 4 . 

We have chosen MnW0 4 as an example since it has Mn 2+ ions in polar surroundings 
and is directly applicable to our approach. At the same time magnetoelectrics with 
magnetic ions in noncentrosymmetric environment are numerous and a brief review of 
the recently discovered multiferroics gives the following examples. Similar to MnW0 4 
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the local C2 symmetry is found for Fe 3+ ions in NaFeSi20g [20] and for Ni 3+ "spine" 
spins in N13V2O8 [21] . Local C s symmetry is found for Cu 2+ ions in LiCu202 [22J and 
for one of the Cr 3+ positions in a-CaCr204 [23]. All of the nonequivalent Fe 3+ ions in 
FeVC>4 [21] and FeTe20sBr [25] have local Ci symmetry. For all these magnetic ions 
the linear in r*part of the crystal field expansion can be written as Vcf{t) = Vq + c- f 
with |5| taking values from 9 ■ 10" 10 N for Ni 3 V 2 8 to 6 ■ 1(T 9 N for LiCu 2 2 . 

Magnetoelectric effect was recently found in manganese phosphorus trisulfide 
MnPSs [26J. MnPSs possesses a monoclinic crystal structure C2/m and shows an 
antiferromagnetic collinear order with k = below T^=78 K. Similar to MnWO-4 it has 
two Mn 2+ ions in the unit cell with the local C2 symmetry and our analysis performed 
above for k = magnetic structures in wolframite is directly applicable to MnPSs. 
According to the neutron diffraction data the magnetic structure is characterized by A x 
and A z 7^ below TV- Thus, following the Eqs. (CQ) one can expect linear magnetoelectric 
effect below T N with magnetoelectric interactions P^A x F y , P^A z F y , P y A^F^, P y A x F z 
and P y A z F x , where \l — x,z. According to Eqs. ([6]) the magnetoelectric coefficient 
a yz = dP y /dH z , for example, can be estimated as a yz = d 0z (l/v )(X/A) 2 A Z ■ dF z /dH z , 
with dF z /dH z ~ 2 • 10~ 6 Oe" 1 [27] and A z ~ 1. The polar local distortion of the Mn 2+ 
environment in MnPS3 is much smaller than in wolframite giving c z ~ 1.4 • 10 -11 N. 
Using v ~ 207 A 3 we obtain a yz ~ 3.4 • 10~ 8 /iC-m^-Oe -1 , which is rather small. 

LiNiP04 possesses an orthorhombic symmetry with space group Prima and 
shows linear magnetoelectric effect in the low temperature C-type commensurate 
antiferromagnetic phase with k = below 20.8 K [281 EH] - The spins are predominately 
directed along the c axis and the magnetic structure is described by the order parameter 
C z transforming according to the IR GM 4_ . The phenomenological magnetoelectric 
interactions C Z M Z P X and C Z M X P Z were suggested earlier [28] and here we can estimate 
the magnetoelectric coefficient according to our microscopic model. Our microscopic 
approach differs from that suggested for LiNiPC>4 earlier [30], which is based on 
lowering the superexchange interaction energy due to the uniform displacement of 
oxygen tetrahedra. The local symmetry of the Ni 2+ ions surroundings is C s with 
|5| » 4.3- 10" 9 N. Similar to the above case of MnPS 3 using dM x /dH x ~ 2-10~ 2 /i B /T per 
Ni-atom [29] we obtain an estimation of the value of magnetoelectric coefficient a zx ~ 
0.31 /iC-m _2 -T _1 in good agreement with the experimental value of 0.2 /xC-m _2 -T _1 [28] 
and about three orders of magnitude higher than that in MnPS3. 

At the same time in many other magnetoelectrics magnetic ions are located in 
centrosymmetric positions. Such is the case, for example, in the rare-earth manganites 
RMnOs [31] and CuO [32] where the local symmetry around Mn 3+ and Cu 2+ is Q. 
Nevertheless, our approach is valid also in these cases. The application of our model 
becomes more complicated and will be published elsewhere [33], but briefly can be 
described as follows. The local symmetry around the rare-earth ions in RMnC>3 is 
C s . The importance of rare-earth ions in formation of electric polarization in RMnC>3 
was recently pointed out [311 EH]- Indeed, the magnetic order of Mn 3+ ions induces 
magnetic ordering of the rare-earths through various exchange mechanisms, which makes 
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our approach applicable. In the case of CuO the oxygens have local symmetry C2. 
The role of oxygen in the superexchange is still a subject of debate jTU [36] . Being 
the intermediate ion conducting superexchange, O 2- should possess induced magnetic 
moment when copper spins order, which again allows application of our model. The 
antiferromagnetic spin polarization at the oxygen sites was measured, for example, in 
the multiferroic TbMn 2 05 |37J. 

5. Conclusions 

We have suggested a microscopic model of magnetoelectric interactions, which 
directly exploits the fact that in many magnetoelectrics magnetic ions are located in 
noncentrosymmetric positions. The model is illustrated by the examples of M11WO4 
and LiNiP04, for which we obtained good correspondence of the values of electric 
polarization and magnetoelectric coefficient, respectively. We also give an estimate 
of the magnetoelectric coefficient in MnPS3. 
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